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The balance function is a new observable based on the principle that charge is locally conserved 
when particles are pair produced. Balance functions have been measured for charged particle pairs 
and identified charged pion pairs in Au+Au collisions at ^/snn = 130 GeV at the Relativistic 
Heavy Ion Collider using STAR. Balance functions for peripheral collisions have widths consistent 
with model predictions based on a superposition of nucleon-nucleon scattering. Widths in central 
collisions are smaller, consistent with trends predicted by models incorporating late hadronization. 

PACS numbers: 25.75.Gz 



Collisions of Au nuclei at ultra-relativistic energies can 
produce large energy densities and high temperatures. 
At these densities and temperatures, the produced mat- 
ter may be best portrayed by partonic (quark/gluon) de- 
grees of freedom as opposed to those characterizing a hot 
hadronic phase. This partonic phase would necessarily 
be followed by a transition to normal hadronic particles, 
which are ultimately measured [lj. Numerousprobes of 
the high density medium have been proposed Q , includ- 
ingobservables related to fluctuations and correlations 

QDQDQ 

One such observable, the balance function, is sensi- 
tive to whether the transition to a hadronic phase was 
delayed, as expected if the quark-gluon phase were to 
persist for a substantial time || . The authors of Ref. Q 



formulate the balance function as follows in this para- 
graph. As a result of local charge conservation, when 
particles and their anti-particles are pair produced, they 
are correlated initially in coordinate space. If hadroniza- 
tion occurs early, the members of a charge/anti-charge 
pair would be expected to separate in rapidity due to 
expansion and re-scattering in the strongly interacting 
medium. Alternatively, delayed hadronization would 
lead to a stronger correlation in rapidity between the 
particles of charge/anti-charge pairs in the final state. 
Measuring this correlation involves subtracting uncorre- 
cted charge/anti-charge pairs on an event-by-event ba- 
sis. The remaining charge/anti-charge particle pairs are 
examined to determine the correlation as a function of 
the relative rapidity, Ay, between the members of the 
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pairs. The balance function is expected to be narrower 
for a scenario with delayed hadronization, and is there- 
fore sensitive to the conjecture that a quark-gluon plasma 
may be produced. 

The balance function is a new observable for heavy 
ion collisions. This function is similar in form to the 
observable used in Q to study the associated charged 
density balance from both p+p and e + + e~ collisions 

In this analysis the balance function is used to examine 
the pseudorapidity correlation of non-identified charged 
particles and the rapidity correlation of identified charged 
pions. The balance function is calculated as: 



2 1 A, 
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where A^ in the case of identified charged pion pairs 

denotes the number of identified charged pion pairs in 
a given rapidity range Ay = \y{^ + ) — y(7r~)|, similarly 

for A ++ , A__, and A |_. In the case of non- identified 

charged particle pairs, pseudorapidity (77) is used. The 

terms A_| , A ++ , A__, and A |_ are calculated using 

pairs from a given event and the resulting distributions 
are summed over all events. Specifically the distribu- 
tion A H is calculated by taking in turn each positive 

particle in an event and incrementing a histogram in Ay 
(Ari) with respect to all the negative particles in that 

event. The distribution A_| is then summed over all 

events. A similar procedure is followed for A ++ , A__, 
and A |_. is the number of positive(negative) pi- 

ons or non-identified charged particles summed over all 
events. The balance function is calculated for all charged 
pion or charged particle pairs for a given centrality bin. 
The balance function, in addition to being proportional 
to opposite charge correlations, is also proportional to the 
acceptance for each Ay or A77 bin. Thus the measured 
balance functions presented in this paper are narrower 
than model calculations assuming no detector acceptance 
effects. 

The data used in this analysis were measured using 
the Solenoidal Tracker at RHIC (STAR) detector for 
Au+Au collisions at ^/snn — 130 GeV. The main de- 
tector used for this analysis was the Time Projection 
Chamber (TPC) UK located in a solenoidal magnetic 
field of 0.25 T The TPC provided tracking informa- 
tion for charged particles having transverse momentum 
Pt > 100 MeV/c and \rj\ < 1.3 with complete azimuthal 
acceptance. Minimum bias events were defined by the 
coincidence of two Zero Degree Calorimeters (ZDCs) [3 
located ± 18 m from the center of the interaction region. 
Central collisions were defined using both the ZDCs and 
the Central Trigger Barrel (CTB), an array of scintil- 
lator slats surrounding the outside of the TPC. A co- 
incidence in the ZDCs accompanied by a signal in the 
CTB above a threshold empirically determined to corre- 
spond to approximately the 15% most central events was 



used to define the central collision trigger. Centrality 
was determined using the charged particle multiplicity 
distributions. Central collisions correspond to - 10%, 
mid-central 10 - 40%, mid-peripheral 40 - 70%, and pe- 
ripheral 70 - 96% of the integrated total. Only events 
with five or more tracks in the TPC were used in this 
analysis. 

For this analysis 150,777 minimum bias and 111,142 
centrally triggered events were used. Only charged par- 
ticle tracks having more than 15 space points along the 
trajectory were accepted. In addition the reconstructed 
trajectory was required to point within 1 cm of the pri- 
mary vertex. To suppress double counting due to track 
splitting, the ratio of reconstructed space points to pos- 
sible space points along the track was required to be 
greater than 0.52. Tracks were further required to have a 
momentum 0.1 < p < 2 GeV/c. Simulations of TPC per- 
formance indicate the tracking efficiency for tracks within 
this acceptance and momentum range is 85-90% [lij . 

Charged pion identification was accomplished by se- 
lecting particles within 2 standard deviations of the ex- 
pected dE/dx for a given momentum, which provided 
pion identification for momenta less than 0.7 GeV/c. 
Kaon contamination of the pions is negligible at low mo- 
mentum and is estimated to be 5% at 0.7 GeV/c. Elec- 
trons were removed by requiring the measured dE/dx to 
be more than 2 standard deviations away from the ex- 
pected value for an electron of the measured momentum. 
An estimated 1% of the electron contamination remains 
in the pion sample after these cuts are made. The lowest 
A77 (Ay) bin, which was estimated to contain 90% of the 
remaining electron contamination, was not used in the 
calculation of the width of the balance function. 

Theoretical simulations of the balance function for 
Au+Au collisions at ^/s^n = 130 GeV were carried out 
using version 1.36 of HIJING ^5| with the default set- 
tings for impact parameters ranging from to 15 fm. 
A realistic distribution was used for event vertices within 
the interaction region. HIJING events were processed us- 
ing GEANT and the TPC track reconstruction soft- 
ware. The same cuts were applied to the HIJING events 
that were applied to the data. None of the features of 
the balance functions computed from the filtered simula- 
tions show any dependence on centrality 17]. Thus the 
STAR results are compared to HIJING integrated over 
centrality. Figure^ shows the balance function (Eq. 1) 
measured for charged particles for the central and pe- 
ripheral Au+Au collision samples. The errors shown are 
statistical. The systematic error in the balance function 
is estimated to be 5% due to systematic variations in 
tracking efficiency, the measurement of pseudorapidity, 
and contamination of electrons. The comparable bal- 
ance function derived for HIJING events simulated in 
the STAR detector using GEANT is shown in Figure^D. 

The area under the balance function is constrained by 
global charge conservation and STAR's acceptance |18| . 
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FIG. 1: The balance function versus Ar/ for charged parti- 
cle pairs from a) central and peripheral Au+Au collisions at 
-^/sjvjv = 130 GeV and mixed events from central and periph- 
eral Au+Au collisions, and b) HIJING events filtered with 
GEANT and shuffled pseudorapidity events from central 
and peripheral Au+Au collisions. To guide the eye, Gaus- 
sian fits excluding the lowest bin in Ar/ are shown. The error 
bars shown are statistical. The balance function for HIJING 
events is independent of centrality. 



Physical effects over and above this constraint can be 
discerned by comparison to a reference data set that pre- 
serves global charge conservation, while removing effects 
of dynamical particle correlations. A relevant reference is 
provided in Figure lb for central and peripheral collision 
samples independently by calculating the balance func- 
tion after the pseudorapidities of all charged particles 
within each measured event have been randomly shuf- 
fled. Dynamical correlations in Au+Au are reflected in 
the deviation of the results in Figure la from the shuffled 
pseudorapidity results in Figure lb. In addition Figure 
la also shows the balance functions generated from con- 
ventional mixed-event samples constructed by choos- 
ing random particles from different measured events with 
similar event vertex positions and centralities. The bal- 
ance function for mixed events integrates to zero because 
global charge conservation has been removed. The fact 
that the balance function is zero for all Ar/ for all cen- 
tralities demonstrates that STAR's acceptance in A?/ is 
smooth. For both the mixed events and shuffled pseudo- 
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FIG. 2: The width of the balance function for charged par- 
ticles, (Ar/), as a function of normalized impact parameter 
(b/bmax)- Error bars shown are statistical. The width of 
the balance function from HIJING events is shown as a band 
whose height reflects the statistical uncertainty. Also shown 
are the widths from the shuffled pseudorapidity events. 



rapidity samples, the measured inclusive pseudorapidity 
distributions are preserved. 

Within this area constraint, the variation of the bal- 
ance function with centrality can be effectively charac- 
terized by the single parameter (At/), the mean pseudo- 
rapidity difference weighted by the balance function (ex- 
cluding the lowest bin in Ar] to reduce the background 
correlation from electron contamination). We refer to 
(Ar/) below as the "width" of the balance function. The 
measured widths for four centrality classes are shown in 
Figure as a function of the impact parameter fraction 
b/b m ax, which is determined using a simple geometrical 
picture to relate impact parameter to fractions of the 
total cross section. In Figure [2J the width of the balance 
function measured for central collisions is significantly 
smaller than that for peripheral collisions. The results 
for the mid-peripheral and mid-central centrality classes 
decrease smoothly and monotonically from the periph- 
eral collision value. Figure [21 indicates that while the 
width observed in peripheral collisions is consistent with 
the HIJING prediction, the balance function for central 
collisions is significantly narrower, suggesting a variation 
in the underlying particle production dynamics between 
these two classes of events. In Figure 2 the widths from 
the shuffled pseudorapidity events are also shown. These 
widths show little centrality dependence and are wider 
than those of the data or HIJING. The widths from shuf- 
fled pseudorapidity events represent the maximum possi- 
ble width of a balance function measured with the STAR 
detector. 

The results for identified charged pion pairs are similar 
to those for non-identified charged particles as indicated 
in Figure |31 The overall shape of the balance function is 
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FIG. 3: The balance function versus Ay for identified 
pion pairs from a) central and peripheral Au+Au collisions 
at ^/snn = 130 GeV, and b) HIJING events simulated in 
GEANT [l6j|. To guide the eye, Gaussian fits excluding the 
lowest two bins in Ay are shown. The error bars shown are 
statistical. The balance function for HIJING events is inde- 
pendent of centrality. 



similar to that in Figure ^ However, the data for pions 
have a dip near Ay — (Figure . As shown in ref- 
erence [I8j, this dip can be understood as the combined 
effect of Bose-Einstein correlations and Coulomb interac- 
tions between charged pions. HIJING does not account 
for these effects, and the balance function predicted in 
Figure |2b therefore does not show a dip, although the 
enhancement of the lowest bin due to electron contami- 
nation is still apparent. 

The width of the balance functions for the four cen- 
trality classes is shown in Figure 0] where the lowest two 
bins in Ay have been excluded to avoid the effects of 
Bose-Einstein correlations and Coulomb interactions and 
enable a valid comparison with HIJING. The results in- 
dicate that the width of the balance function for periph- 
eral events is consistent with that expected from HIJING, 
while the width for central events is significantly smaller 
as was observed with charged particle pairs. 

To gain some insight into the observed widths in cen- 
tral collisions, we compared with the thermal model (fi- 
nal temperature of 165 MeV) presented in || filtered 
through the acceptance of STAR. We find that the pre- 
dicted widths are larger than those we observe in cen- 
tral collisions. Thermal model calculations have also 
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FIG. 4: The width of the balance function for identified 
charged pions, (Ay), as a function of normalized impact pa- 
rameter (b/bmax)- Error bars shown are statistical. The width 
of the balance function from HIJING events is shown as a 
band whose height reflects the statistical uncertainty. 



been done at low final temperature (105 MeV) and high 
transverse velocity (0.77 c), while maintaining the same 
average transverse momentum. These calculations pre- 
dict a balance function width consistent with our ob- 
servations in central collisions. Thus the observed nar- 
rowing of the balance function in central collisions can 
be parametrized in a thermal model incorporating both 
strong transverse flow and a constraint that the balanc- 
ing particles are emitted close together in space-time [i^ • 
This constraint could arise from delayed hadronization 
compared with the characteristic 1 fm/c hadronization 
time or from some other phenomena such as anomalously 
short diffusion of particle pairs [i(J ■ Other effects besides 
flow may also narrow the balance function with central- 
ity. One effect is the decay of resonances such as the 
p° that have lifetimes similar to the proposed time of 
delayed hadronization. Another effect could be opacity 
|'2l| where pairs with large opening angles are absorbed 
leading to a narrower balance function. 

In summary, measurement of the balance function for 
Au+Au collisions at GeV has been stimu- 

lated by the prediction [8j that the width of the bal- 
ance function should be significantly reduced by late 
hadronization. We indeed observe a narrowing of the bal- 
ance function for more central collisions for all charged 
particle and for charged pion pairs. Only for peripheral 
collisions is the width consistent with HIJING predic- 
tions treating the Au+Au collision as a superposition of 
independent nucleon-nucleon scatterings. Interpretation 
of the observed narrowing requires more detailed study 
of its sensitivity to such other effects as flow, resonance 
production and opacity, in addition to late hadronization. 
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